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Synopsis 


An  important  method  of  investigating  the  physical  properties  of 
high-temperature  laboratory  and  solar  plasmas  is  the  analysis  of 
optically  thin  spectral  lines  emitted  by  multiply- charged  ions.  In 
addition,  the  energy  loss  due  to  the  emitted  far  ultraviolet  and 
x-ray  radiation  may  have  a significant  influence  on  the  dynamical 
behavior  of  the  plasma.  In  order  to  understand  and  interpret  the 
optical  emissions  from  these  plasmas  requires  a knowledge  of  the 
various  processes  contributing  to  the  observed  spectra.  One  such 
process,  and  oftentimes  the  dominant  one,  is  electron  impact 
excitation  of  highly-charged  ions.  Because  of  the  paucity  of  both 
experimental  and  theoretical  cross  section  data  for  multiply-charged 
ions,  there  is  an  overwhelming  demand  for  this  information.  In 
this  paper  we  have  calculated  a number  of  selected  transitions  in 
multicharged  silicon  ions  from  Si  VI  to  Si  XIV.  The  calculations 
were  carried  out  using  the  method  of  distorted  waves  ignoring  the 
effects  due  to  exchange  except  for  the  excitations  in  heliumlike 
Si  XIII,  where  it  was  included.  Included  in  the  calculations  are 
values  for  the  line  strength,  collision  strength  and  excitation 
rates.  In  audition,  some  of  these  rate  have  been  represented  by 
a simple  two  parameter  fit  for  convenience. 

This  paper  is  being  submitted  for  publication  to  the  Journal 
of  Quantitative  Spectroscopy  and  Radiative  transfer. 
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Abstract 

Line  strengths,  collision  strengths  and  excitation  rates  have 
been  calculated  for  a variety  of  transitions  in  multicharged  silicon 
ions  from  Si  VI  to  Si  XIV.  The  collision  strengths  were  evaluated 
in  an  LS  coupling  scheme  in  the  distorted  wave  approximation 
neglecting  exchange  except  for  the  helium-like  transitions. 
Excitation  rates  were  then  obtained  by  integrating  the  collision 
strength  over  a Maxwellian  velocity  distribution  function.  These 
results  are  then  described  by  a simple  two-parameter  fit  for 
the  rates. 
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I • Introduction 

The  radiative  emissions  produced  from  high  temperature  Laboratory  and 
solar  plasmas  contains  a potential  wealth  of  information  on  such  parameters 
as  chemical  composition,  temperature  and  density  structure,  and  a measure 
of  the  radiative  cooling.  However,  before  any  reliable  quantitative  infor- 
mation can  be  deduced  from  the  observed  spectra,  some  knowledge  must  be 
available  on  the  fundamental  processes  leading  to  its  formation.  One  such 
process,  and  often  the  dominant  one  in  an  optically  thin  plasma,  is 
electron -impact  excitation  of  highly  charged  ions.  Since  a large  fraction 
of  the  emitted  line  radiation  is  a direct  consequence  of  the  interaction  of 
electrons  and  multi-charged  ions,  it  is  essential  that  accurate  values  of 
the  cross  sections  and  excitation  rates  be  available  for  the  interpretation 
of  the  observed  spectra.  Because  of  the  paucity  of  both  experimental 
and  theoretical  cross-section  data  for  multi-charged  ions,  there  is  an  over- 
whelming demand  for  this  information. 

In  an  earlier  paper, we  provided  values  of  the  cross  sections  and 
excitation  rates  for  a number  of  transitions  in  various  ionization  stages 
for  several  ions  of  solar  interest:  in  particular,  iron,  neon,  magnesium,  and 
aluminum.  Maintaining  the  spirit  of  this  paper,  we  will  present  similar 
results  for  multi-charged  silicon  ions.  In  an  ideal  situation, it  would  make 
sense  to  obtain  these  data  by  employing  the  best  available  method,  namely,  the 
close-coupling  method  with  correlation.  However,  this  would  be  impractical 
in  the  case  of  many  transitions.  On  the  other  hand,  the  distorted-wave  method 
(DW)  is  appropriate  for  ionized  systems  and  has  an  anticipated  accuracy  that  in- 
creases as  the  charge  on  the  target  ion  increases;  it  is  also  suitable  for  dealing 
with  a number  of  transitions  in  an  economical  fashion,  a consideration  that 
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can  not  be  ignored. 


II.  Method  of  Calculation 


The  method  of  distorted  i&ves  was  used  to  obtain  collision  parameters 
for  all  transitions.  In  an  IS  coupling  scheme,  the  excitation  cross  section 
between  atomic  terms  a LS  -*  Gf*L  'S  * is  written  as 

Q (aLS  -*  a*L 's ')  « 2Lh  s ) v a2  , (1 

E(2L  + 1 ) (2S+1  ) ° 

where  E is  the  energy  in  Rydbergs  of  the  incident  electron  and  fj  (aLS,  a'L'S'') 
is  the  collision  strength  which  is  symmetric  with  respect  to  aLS  and  a'L 's  '. 

A discussion  of  the  method  used  in  these  calculations  is  presented  in 
Ref.  (1)  and  we  refer  the  reader  there  for  the  details.  As  in  Ref.  (1), 
we  made  the  following  simplifying  assumptions:  Wave  functions 

of  the  target  ions  were  represented  by  I.S  coupling  and  we  have  omitted  effects 
of  configuration  mixing,  the  electron  exchange  was  taken  into  account  only 
for  collisions  with  helium-like  Si  XIII,  in  the  expansion  of  the  expression 
for  the  interaction  potential,  we  have  retained  only  terms  corresponding  to 
multipole  orders  lower  than  _“>. 

Although  the  collision  strength  is  the  fundamental  atomic  parameter, 
the  result  required  for  the  interpretation  of  observed  spectra  is  the 
collision  strength  averaged  over  an  appropriate  velocity  distribution.  For 
plasma  electrons  described  by  a Maxwellian  velocity  distribution  function 
f(v)  and  electron  temperature  T,  the  excitation  rate  coefficient  can  be 
written  as 


xae  «> 


-I 


f (V) 


v Qa8  (V) 


cm  V sec . 


(2) 


In  terms  of  the  collision  strength,  the  rate  becomes 


c 00 

Xap  (T)  = f Q exp  (-E/kT)  d (E/kT) , 

h*  T*  ~ 


Eo/kT 


where  and  E is  the  statistical  weight  of  level  ar  and  E is  the 
(X  o o 


excitation  energy.  Results  will  be  presented  for  both  the  collision 


strength  and  excitation  coefficient. 


III.  Results 


The  charge  states,  transitions,  wavelengths  and  line  strengths 


for  which  we  have  performed  calculations  are  presented  in  Table  1.  The 


target -ion  ground -state  wave  functions  are  of  the  Hartree-Fock  type,  as 


given  by  Clementi.  The  excited -state  wave  functions  are  obtained 


from  a semi-empirical  method  using  experimental  eigenvalues  when  avail - 


able;  otherwise,  theoretical  estimates  are  used.  The  method  generates 


radial  wave  functions, which  are  obtained  in  the  field  of  a screened  atomic 


potential  and  a statistical  exchange  potential  with  a variable  parameter 


that  guarantees  the  proper  asymptotic  behavior.  In  order  to  ensure  that 


the  radial  functions  satisfy  the  orthogonality  condition,  the  method  of 


undetermined  Lagrange  multipliers  is  utilized  in  a self-consistent  fashion 


with  the  radial  function  routine. 


After  generating  the  excited -state  wave  functions,  we  formed  the 


dipole  matrix  elements  (or  dipole-transition  integral)  and  obtained 


estimates  for  the  line  strength.  These  results  are  presented  in  Table  1 


along  with  some  values  quoted  by  Wiese  et  al  ' given  in  parentheses. 


For  those  transitions  where  a comparison  could  be  made, the  agreement  is 


j 


reasonably  good.  Differences  are  probably  caused  by  the  effect  of  con- 
figuration mixing, which  was  ignored  in  our  calculation.  In  addition,  we 
have  also  listed  the  quadrupole-transition  integral  where  appropriate. 

For  convenience,  we  give  the  relationship  between  the  spontaneous  decay  rate 
A and  the  line  strength  S (atomic  units)  for  dipole  transitions  only 
[see  Ref.  0)  for  quadrupole  transitions,  etc.],  viz. 

A = -i.0^6.  x 10 — s (sec  ^ 0) 

uJ  A5 

where  u)  is  the  statistical  weight  of  the  upper  level  and  X is  the  wave- 
length in  angstrom  units. 

The  results  of  the  collision  strength  calculations  are  presented  in 
Figs.  1 through  7 as  functions  of  incident  electron  energy  normalized 
with  respect  to  the  excitation  energy,  Eq,  of  the  transition.  The 
curves  in  the  various  figures  are  numbered  in  accordance  with  their  entry 

in  Table  1.  For  example,  curve  1 in  Fig.  1 refers  to  the 

5 2 h ,3  2 

2p  P-2p  ( P)^d  D transition  of  Si,  VI  and  is  identified  as  such  by 
referring  to  the  column  labeled  "CODE"  in  Table  1.  A review  of  the 
figures  indicates  that  the  collision  strength  is  fairly  constant  for  the 
An  = 0 transitions  but  varies  considerably  otherwise,  particularly  for 
the  spin  exchange  excitations. 

Some  caution  must  be  exercised  in  interpreting  the  results  presented 
in  Table  1.  The  designation  of  a J value  in  the  description  of  a tran- 
sition is  used  only  to  identify  the  value  of  the  energy  level  used  in 
the  calculation.  We  have  not  calculated  collision  parameters  for  J-J' 
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transitions.  The  results  in  the  Table  refer  to  the  multiplet  average. 
However,  the  results  can  be  transformed  easily  to  yield  such  information, 
provided  that  we  neglect  small  energy  differences  of  J sublevels. 
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For  dipole  transitions , 


Cl  f SI., si.  *1  <2J  + 1>  (2J'  + l)  , (5) 

2S  + 1 


_ , , j 1 J J'  I 

n (SLJ.SL  J ) = j sl'l  j 


and,  if  S = 0 or  L = 0,  then  generally 


('SLJ.S'l'j')  = __l£_ Li n (SL.S'L'). 

(2S  + 1)  (2L  + 1) 


The  excitation  rates  were  obtained  from  eqn.  (j) . For  regions 
exceeding  the  range  of  our  calculation,  Q was  assumed  to  be  constant. 
This  procedure  introduces  only  a small  error  into  the  excitation  rate. 

The  majority  of  rate  calculations  can  be  represented  adequately  by 
a parametric  fit  of  the  form 


X = a e 


-0/T 


n-V 


(6) 


where  (7,  8 and  V are  the  fitted  parameters  and  the  electron  temperature  T is 
in  ev.  We  have  divided  the  results  into  two  groups.  In  the  first  group,  we  have 
set  V = with  the  values  for  a.  and  8 shown  in  Table  1.  The  fit  was  con- 
sidered good  when  it  reproduced  the  calculated  values  to  within  20$  over 
a temperature  range  oflOOtolc  kev.  The  second  group  of  rates  are  pre- 
dominantly intersystem  transitions.  Since  a single  value  of  V would  not 
represent  the  rate  over  the  temperature  range  of  interest,  the  results  are 
presented  graphically  in  Figs.  M and  0. 


IV.  Summary 


We  have  calculated  line  strengths,  collision  strengths  and  rate 
coefficients  for  a variety  of  transitions  in  multicharged  silicon  ions 
from  Si  VI  to  Si  XIV.  The  line  strengths  were  obtained  by  using  dementi 
wave  functions  for  the  ground-state  configurations  and  excited -state 
wave  functions  generated  by  a semi-empirical  method.  For  those  cases 
where  comparisons  could  be  made, the  agreement  was  good.  The  collision 
strengths  were  calculated  in  an  LS  coupling  scheme  in  the  distorted- 
wave  approximation,  neglecting  exchange  except  for  the  helium-like 
transitions.  These  results  were  then  integrated  over  a Maxwellian 
velocity  distribution  function  to  yield  rate  coefficients.  The  rates 
were  presented  either  graphically  or  in  terms  of  a two- parameter  fit. 
Because  of  the  paucity  of  both  experimental  and  theoretical  data  for 
multicharged  ions,  we  have  been  unable  to  make  comparisons  with  other 
studies  to  judge  the  accuracy  of  our  calculations.  These  results  will 
be  compared  in  a forthcoming  paper  with  those  currently  under  investigation 
using  multiconf igurational  wave  functions. 


J.  Davis,  P.  C.  Kepple,  and  M.  Blaha,  JQSRT  16,  10^.5  (1077) 
E.  Cleraenti,  IBM  Journal  of  Res.  and  Dev.  Q 2 (1065) 


3.  W.  Wiese,  M.  Smith, and  B.  Miles 


Atomic  Transition  Probabilities  " 


NBS  22,  Vol . 1 (I966) 


•m 


Appendix 


In  a previous  paper^  the  following  misprints  should  be  corrected: 


Table  3,  the  last  line:  the  values  8.9  and  1230.0  should 
be  replaced  by  8.O(-10  and  12300.0,  respectively. 


Equation  (21),  the  last  line:  l in  the  second  3-j  symbol 

.1 


should  be  replaced  by  A 


The  titles  of  Figs.  7 and  8 should  be  interchanged. 


Equation  (18)  symbol  L should  be  replaced  by  1^  , 

A.  \ * 
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